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Abstract 

Assume that samples of a filtered version of a function in a shift-invariant space 
are avalaible. This work deals with the existence of a sampling formula involving 
these samples and having reconstruction functions with compact support. Thus, low 
computational complexity is involved and truncation errors are avoided. This is done 
in the light of the generalized sampling theory by using the oversampling technique: 
more samples than strictly necessary are used. For a suitable choice of the sampling 
period, a necessary and sufficient condition is given in terms of the Kronecker canonical 
form of a matrix pencil. Comparing with other characterizations in the mathematical 
literature, the given here has an important advantage: it can be reliable computed 
by using the GUPTRI form of the matrix pencil. Finally, a practical method for 
computing the compactly supported reconstruction functions is given for the important 
case where the oversampling rate is minimum. 
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1 Statement of the problem 

Let be a shift-invariant space in L 2 (M) with stable generator (p € L 2 (R), i.e., 
V v := {f{t) = J>n <p(t - n) : {a n } G £ 2 (Z)} C L 2 (R) . 
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Nowadays, sampling theory in shift-invariant spaces is a very active research topic (see, for 
instance, [1, 2, 3, 4, 5, 10, 21, 22] and the references therein) since an appropriate choice 
for the generator ip (for instance, a B-spline) eliminates most of the problems associated 
with the classical Shannon's sampling theory [18]. 

Suppose that a linear time- invariant system C is defined on V v . Under suitable condi- 
tions, Unser and Aldroubi [4, 17] have found sampling formulas allowing the recovering of 
any function / G V v from the sequence of samples {(£/) (n)} ne z. Concretely, they proved 
that for any / G V^, 

f(t) = J2^f(n)S c (t-n), t€R, (1) 

where the sequence {Sc(t — n)} ne z is a Riesz basis for V^. Note that a reconstruction 
function Sc with compact support implies low computational complexity and avoids trun- 
cation errors. Even when the generator <p has compact support, rarely the same occurs 
with the reconstruction function Sc in formula (1). A way to overcome this difficulty is 
to use the oversampling technique, i.e, to take samples {(£/) (wT)}„ £ 2 with a sampling 
period T < 1. This is the main goal in this paper: Assuming that both the generator p 
and C<p have compact support, we study the existence of stable sampling formulas with 
compactly supported reconstruction functions, which allow us to recover any / G V v from 
the samples {(£/) (nT)} n6 z, where the sampling period is T := r/s < 1 for some positive 
integers r < s. This is done in the light of the generalized sampling theory obtained in 
[12] by following an idea of Djokovic and Vaidyanathan in [10]. For the sake of notational 
simplicity we have assumed that only samples from one linear time-invariant system £ are 
avalaible. 

In so doing, the problem is connected with the search of polynomial left inverses of a 
certain s x r polynomial matrix G{z) intimately related to the sampling problem. Taking 
advantage of the special structure of the matrix G(z) we give a necessary and sufficient 
condition which involves the matrix pencils theory. Concretely, this condition uses some 
information contained in the Kronecker canonical form of a matrix pencil associated with 
the matrix G(z). From a practical point of view, this information can be retrieved from the 
GUPTRI (General UPper TRIangular) form of the matrix pencil. It is worth to mention 
that the GUPTRI form can be stably computed. 

The mathematical problem of finding a polynomial left inverse of a polynomial matrix 
G(z) has been studied in [7] by Cvetkovic and Vetterli in the filter banks setting. It 
involves the Smith canonical form S(z) of the matrix G(z). Roughly speaking, the Smith 
canonical form of G(z) must contain monomials in its diagonal. From a practical point of 
view, the Smith canonical form has an important drawback: there is not a stable method 
for its computation. 

Another algebraic approach is the following (see, for instance, [16]): Assume that G(z) 
is a s x r Laurent polynomial matrix (r < s); whenever the greatest common divisor of 
all minors of maximum order r is a monomial, then its Smith canonical form S(z) has 
monomials in its diagonal. Without loss of generality we can assume that the 7 := (^) 
minors of order r in G(z) are polynomials with positive powers in z. Invoking Euclides 
algorithm we can obtain (*) polynomials, fi(z), . . . , fj(z), such that 

7 

fn(z)A n (z) = m(z) , for all z G C 

71=1 
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where A n , 1 < n < 7, are the minors of order r of G(z) and m(z) is a monomial. Denote by 
D' n (z) the adjoint matrix corresponding to the minor A n and D n (z) the matrix obtained 
from D' n (z) by adding s — r zero columns. Thus, Z) n (z)G(z) = A n (z)l r , and consequently 

(f^(z)A^)W) = I r , 

where /^(^) := f(z)/m(z) could be a Laurent polynomial, 1 < n < 7. From a practical 
point of view the drawback here is the effective calculation of the (*) minors of G(z) 
whenever r becomes larger. 

The paper is organized as follows: In Section 2 we include the needed preliminaries to 
understand the raised problem. The existence of reconstruction functions with compact 
support depends on the rank, for z G C \ {0}, of a polynomial matrix G(z), associated 
with the sampling problem. In Section 3, a suitable choice of the sampling period T = r/s 
reduces our problem to a matrix pencil problem. Thus, we give a necessary and sufficient 
condition for the existence of compactly supported reconstruction functions which involves 
the Kronecker canonical form of a singular matrix pencil. Section 4 is devoted to compute 
a polynomial left inverse of the matrix G(z) in the important case where the oversampling 
rate is minimum. Finally, we briefly remind, as an Appendix, the canonical forms alluded 
in what follows. 

2 Preliminaries on generalized sampling 

From now on, the function tp G L 2 (R) is a stable generator for the shift-invariant space 
V v ■= {/(*) = an ~ n ) : M G £ 2 (%)} C L 2 (R) , 

i.e., the sequence {ip(- — n)} n< =z is a Riesz basis for V v . A Riesz basis in a separable Hilbert 
space is the image of an orthonormal basis by means of a bounded invertible operator. 
Recall that the sequence {ip(- — n)} n& z is a Riesz basis for if and only if 

< ||$||o < Halloo < OO, 

where ||$||o denotes the essential infimum of the function $>(w) := J2kez \<P( W + ^)| 2 i n 
(0, 1), and H^Hoo its essential supremum (ip denotes, as usual, the Fourier transform of tp). 
Furthermore, \\&\\o and H^Hoo ar e the optimal Riesz bounds [6, p. 143]. 

We assume throughout the paper that the functions in the shift-invariant space V v 
are continuous on R. Equivalently, the generator <p is continuous on R and the function 
Snez Iffy ~ n )\ 2 i s uniformly bounded on R (see [22]). Thus, any / G is defined as 
the pointwise sum f(t) = ^ n ^iO, n (p{t — n) on R. Besides, is a reproducing kernel 
Hilbert space where convergence in the L 2 (R)-norm implies pointwise convergence which 
is uniform on R (see [12]). 

The space V v is the image of L 2 (0, 1) by means of the isomorphism : L 2 (0, 1) — > 
which maps the orthonormal basis {e~ 27rmw } n ^z for L 2 (0, 1) onto the Riesz basis {<p(t — 
n)} n £z for Vtp. Namely, for each F G L 2 (0, 1) the function T^F G V v is given by 

(T v F)(t):=^(F(.),e- 2 ™-) L2(0)1) ^-n), teR. (2) 
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Suppose that C is a linear time-invariant system denned on V v of one of the following 
types (or a linear combination of both): 

(a) The impulse response h of C belongs to n L 2 (R). Thus, for any / G V v we have 

/oo 
f(x)b(t-x)dx, tel. 
-oo 

(b) C involves samples of the function itself, i.e., (Cf)(t) = fit + d), t G R, for some 
constant d G R. 

For fixed positive integers s > r, consider the sampling period T := r/s < 1. The 
first goal is to recover any function / G Vu> by using a frame expansion involving the 
samples { (£/)(rn/s)} ftgZ . This can be done in the light of the generalized sampling 
theory developed in [12]. Indeed, since the sampling points rn/s, n G Z, can be expressed 
as 

{ rn M„ e z = i rm + C? " 1 ) r / s } me z, i =i,2,...,s ' 
the initial problem is equivalent to the recovery of / G from the sequences of samples 

{Cjf(rn)} ne z,j=i,2,..., s 

where the linear time-invariant systems Cj, j = 1, 2, . . . , s, are defined by 

(£,-/)(*) :=(£/) [i+(i-l)r/ S ], tel. 

Following the notation introduced in [12], consider the functions <7j G L 2 (0, 1), j = 
1, 2, . . . , s, defined as 



9j( w ) ■= ^2(£<P)[n+ (j - l)r/s]e 



—2-rrinw 



(3) 



the s x r matrix 



GH := 



5i (u>) + 
g 2 (w) g 2 {w + \) 



_g s (w) g s {w + \\ 
and its related constants 



i — l ~ 

r - 



j = l,2,...,8 

fc=l,2,...,r 



«g := essinf A min [G*(u;)G(u;)], /3q := ess sup A max [G*(u;)G(t(;)] , 

i«e(0,l/r) u>6(0,l/r) 

where G*(w) denotes the transpose conjugate of the matrix G(w), and A m j n (respectively 
Amax) the smallest (respectively the largest) eigenvalue of the positive semidefinite matrix 
G*(w)G(w). Notice that in the definition of the matrix G(w) we are considering the 
1-periodic extensions of the involved functions gj, j = 1, 2, . . . , s. 

Thus, the generalized sampling theory in [12] (see Theorem 1, Theorem 2 and its proof) 
gives the following sampling result in V v : 
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Theorem 1 Assume that the functions gj, j = 1, 2, . . . , s, defined in (3) belong to L°°(0, 1) 
(this is equivalent to (3q < oo). Then the following statements are equivalent: 

(i) a G > 

(ii) There exist functions a,j in L°°(0, 1), j = 1, 2, . . . , s, such that 

[ ai (w),...,a s (w)] G(w) = [1,0,..., 0] a.e. in (0, 1) . (4) 

(Hi) There exists a frame for V v having the form {Sj(- — rn)} n ^z,j=i,2,...,s such that, for 
any f £ V^, we have 

s 

/ = EE( £ /)[ rn + ^- 1 ) r / s ]^('-™) mL 2 (R). (5) 

neZ j=l 

In case the equivalent conditions are satisfied, the reconstruction functions in (5) are given 
by Sj = rT^aj, j = 1,2, ... ,s, where the functions aj, j = 1,2, ... ,s, satisfy (4). The 
convergence of the series in (5) is also absolute and uniform on R. 

Recall that a sequence {//J is a frame for a separable Hilbert space H if there exist 
two constants A, B > (frame bounds) such that 

Aimi 2 ^!^)! 2 ^!!/!! 2 for ail fen. 
k 

Given a frame {fk} for H the representation property of any vector / G H as a series 
/ = J2k °kfk is retained, but, unlike the case of Riesz bases, the uniqueness of this rep- 
resentation (for overcomplete frames) is sacrificed. Suitable frame coefficients which 
depend continuously and linearly on / are obtained by using the dual frames {gk} of {fk}, 
i.e., {g k } is another frame for H such that / = J2 k (f, gk)fk = fk)9k for each f eH. 

For more details on the frame theory see the superb monograph [6] and the references 
therein. 

It is worth to mention that whenever the 1-periodic functions gj, j = 1, 2, ... , s, are 
continuous on M, the conditions in Theorem 1 are also equivalent to the condition recently 
introduced in [13, Corollary 1]: 



(iv) rankG(tu) = r for all w £ M . 



The goal in this paper is to obtain necessary and sufficient conditions assuring that 
we can find reconstruction functions Sj, j = l,2,...,s, in formula (5) having compact 
support. To this end, assume from now on that the generator (p and dp are compactly 
supported. We introduce the s x r matrix 



G(z) :- 



gl (z) gl (Wz) 

g2(z) g2(Wz) 

.&(*) g s {Wz) 



giiW^z 

g2(W 



r-1. 



gsiW^Z 



(6) 
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where W := e~ 2ni/r and gj(z) := Ym&i&P) [ n + U ~ W 5 ] 3 = 1. 2 ■ ■ ■ , s - Note that 
the matrix G(z) has Laurent polynomials entries, and G(w) = G(e~ 2mw ). On the other 
hand, if the functions 3j(z), j = 1, 2 . . . , s, are Laurent polynomials satisfying 



[ai(z),...,a a (z)]G(z) = [1,0,..., 0] 



(7) 



then, the trigonometric polynomials a,j(w) = 3j(e~ w ), j = 1,2, ...,s, satisfy (4). In 
this case, the corresponding reconstruction functions Sj, j = 1,2, ...,s, have compact 
support. Indeed, in terms of the coefficients Cj >n of 3j(z), that is, aj(z) = Yln& c j,nZ n , 
j = 1, 2, . . . , s, the reconstruction function Sj, j = 1, 2, . . . , s, can be written as (see (2)): 



s j(t) = r '^2 c j,n l P(t-n), te 



(8) 



The existence of polynomial solutions of (7) is equivalent to the existence of a left 
inverse of the matrix G(z) whose entries are polynomials. This problem has been studied 
in [7] by Cvetkovic and Vetterli in the filter banks setting. By using the Smith canonical 
form S(z) of the matrix G(z) (see Appendix 6.1), a characterization for the existence of 
polynomial solutions of (7) has been found in [14]. Namely, assuming that the generator 
if and C(p have compact support, there exists a polynomial vector [ai(z), a2(z), • ■ ■ , a s (z)] 
satisfying (7) if and only if the polynomials ij(z), j = 1,2, ... ,r, on the diagonal of the 
Smith canonical form S(z) of the matrix G(z) are monomials. Assume that the s x r 
matrix 

'h(z) ••• 
i 2 (z) ••• 



S(z) 

















i r (z) 








(9) 



is the Smith canonical form of the matrix G{z) (note that it is the case whenever qg > 0) 
and consider the unimodular matrices V(z) and W(z), of dimension s x s and r x r 
respectively, such that G(z) = V(z)S(z)W(z). 

Observe that if S(z) is the Smith form of the matrix G(z) then, taking into account 
that V(z) and W(z) are unimodular matrices, we have 

rankS(z) = rankG(z) for all zGC. 

Therefore, it is straightforward to deduce that, for each j = 1,2, ...,r, the polynomial 
ij(z) is a monomial if and only if rankS(z) = r for all z G C \ {0}. This condition, under 
the above hypotheses on ip and dp, is equivalent to saying that 



rankG(^) = r for all z G C \ {0} . 



(10) 



(See [14] for the details). The main aim in this work is to search for an equivalent condition 
to (10) useful from a practical point of view. 
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3 Searching for an useful equivalent condition 

The first step is to reduce our polynomial matrix G(z) to a matrix pencil in order to 
use the well-established theory. In so doing we need some preliminaries. Let f(z) = 

a m z m + a m -\z m " 1 H h a\z l + ao be an algebraic polynomial of order m, and let n be a 

positive integer. For each j = 0,1,... ,n—l let fj(z) denote the sum of the monomials a r z r 
where r = j( mod n). Obviously, f(z) = Y^=o fj( z )- The polynomial fj, < j < n — 1, 
is the so-called n-harmonic of order j of the polynomial /; it satisfies: 

g( e 2*»/n z ) = e^^fjiz) for all z e C. 

Assume that supp&p is contained in an interval [0,iV], where JV G N. Thus, the 
functions gj(z) are Laurent polynomials in the variable z. Consider 

r 

p : = minjq G N : q - > 1} . 

s 

It is easy to check that p = c + 1 where c denotes the quotient in the euclidean division 
s\r. Hence, we can write the Laurent polynomials gi(z), j = 1, 2 . . . , s, as: 

gl (z) = £<p(l)z + Cip(2)z 2 + ■■■ + Cip(N - l)z N - 1 

g 2 (z) = &p(^) + C<p(l + ^)z + --- + £<p(N - 1 + ^)z N ~ 1 

g p (z) = <C<p((p ~l) r -)+ C<p{l + (p-l) r -)z + --- + £ip{N-l + (p- l)- s )z N - 1 (11) 
gEp+i(z) = C^{p r - - l)^ 1 + • • • + C^{N - 2+p- s )z N ~ 2 

gs(z) = £<p{(s -l) r --r + l)z-^ + ■■■ + £<p(N - r + 2 + (s - l)^)z N ~ r+2 . 

The polynomial gi(z) has at most N — 1 nonzero terms; the rest of polynomials gj(z), 
2 < j < s, have at most N nonzero terms. In what follows, we use the new matrix 
G(z) = G(z)U(z), where 

U(z) = diag [z^ r - l \ {Wz) { - r - 1 \{W 2 z)^- 1 \. . . , (W- 1 z)( r - 1 )] . 

Thus, all the entries of the polynomial matrix G(z) are algebraic polynomials in z and, 
moreover we have rankG(z) = rankG(z) for all z G C \ {0}. We denote by gj(z) the 
algebraic polynomial z r ~ 1 gj(z), 1 < j < s. 

The strategy is to reduce the polynomial matrix G(z) into another simpler one having 
the same rank for all z G C \ {0}. 

Lemma 1 Consider the matrix G(z) = [Go(z) 62(2) . . . G( r _i)(z)] ; where Gj(z), < 
j < (r — 1), denotes the column vector consisting of the r -harmonics of order j of the 
polynomials gi(z) where 1 < i < s. Then 

G{z) = G{z)n r 

where £l r denotes the Fourier matrix of order r. 
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Proof: For each i = 1, 2, . . . , s we have that gi(z) = J2 T j=o &ij( z ) where gij(z) denotes the 
r-harmonic of order j of gj. We can write the matrix G(z) as 

G(z) = [G (z) + Gi(z) + ■ ■ ■ + G r _i(z) 

G (z) + WG^z) + ■■■ + W^Gr-^z) 



G (z) + W^G^z) + ■■■ + WW 

Hence, in matrix form we have 

G(z) = [Gq(z) Gi(z) . . . G r -i(z)] Q r = G{z)Q r 

where 



1 1 1 

1 W W 2 
1 w 2 W 4 



1 w r ~ l w 2 ^ r ~^ ■ ■ ■ 

is the Fourier matrix of order r. 

Observe that rankG(z) = rankG(z) for all z G C \ {0}. 



1 

W 2( - r ~^ 



□ 



In what follows, we assume that supp&p C [0, iV] and, in addition, we also assume 
that N < r. In this case, having in mind the number of nonzero consecutive terms of the 
polynomial gj(z), we conclude that the r-harmonic of order q, q = 0, 1 ... ,r — 1, of the 
polynomial gi(z), 1 < i < s, is a monomial having the form Ci p z kr+q where Ci q G C and 
k G {0, 1}. This choice of r and, consequently, of the sampling periods T = r/s, r, s <G N 
and s > r, simplifies the structure of the matrix G(z). 

First, let us to give an illustrative example: Consider N = 3, r = 4 and s = 5; here 
T = 4/5,p = 2 and the polynomials gj{z) read 



★z 4 + *z 5 
kz 2 + *z 3 + *z 4 



gl(*) 

g5(-Z) = * + *Z + *Z 2 

Hence, the matrix G(z) reads 



G(z) 



12(2) 

84 («) 



*z 3 + *z 4 + kz 5 
-kz + *z 2 + kz 3 



kZ 4 


kZ 5 








kz 4 


kZ 5 





« 3 


kZ 4 





*z 2 







kZ 


kZ 2 


kZ 3 


k 


kZ 


kZ 2 






(12) 



This example shows that the 3rd and 4th columns have the form z 2 C and z 3 C where 
C,C G C sxl . The first and second columns do not share this property. If we right 
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multiply the matrix G(z) by diag[l, z 1 , 



], we get the new matrix 



G{z) := 



-kZ 
*Z 4 

*z A 




*z 5 


*z 

-kZ 







*z 2 



-kZ 2 





v-2 









* 



0' 

* 

★ * 

* 



Now we can go into the general case for the matrix G(z). Having in mind Eqs.(ll) 
and that gi{z) = z r ~ 1 gi(z) we obtain: 

max { gradgj : 1 < j < s} = (N - 1) + (r - 1) = N + r - 2 < 2r . 
Hence, the matrix G(z) has the form 



G(z) 



'c u z klir c 12 z kl2r+1 ■■■ ci r z kl - r+{r ~^' 



c sl z ksir c s2 z ks2r+l 



k sr r+(r-l) 



where the coefficients kij G {0, 1}. We can easily obtain the following result: 

Lemma 2 Assume that N > 1. Then, for each 1 < j < N — 1 there exist indices i' ^ i, 
1 < i, i' < s, such that kij ^ k^j. Otherwise, for each N < j < r it holds that kij 
for all 1 < i, i' < s. 



ki'j 



Assume that N > 1 and recall that N < r. The entries of the jth column of the matrix 
G(z), where N < j < r, have the form -kz^" 1 (* G C); they could have the form -kz^ 1 or 
^ z r+{j~i) whenever 1 < j < N — 1. Dividing the jth column by z- 7-1 , obviously we obtain 
a matrix with the same rank than G(z) for any z G C \ {0}. Thus, we introduce the new 
polynomial matrix G(z): 

G(z) :=G(z)Q(z) = [M(z)Q] , 

where Q G C sx ( r - N + 1 ) denotes an scalar matrix and Q(z) := diagfl, z^ 1 , . . . , z 1 ~ r ]. When- 
ever vankQ < r — N + 1, we have that rankG(z) = rankG(z) < r for all z G C \ {0} and, 
hence, there is no a polynomial left inverse for G(z). In the case rank^(z) = r — N + 1, 
there exists an invertible matrix R G C sxs such that 

G r 




R 9 = 

where Q' G c^-^ 1 )*^"^ 1 ) is invertible. Thus, 

R G(z) = [RM(z) RQ] -- 



h(z) Q r 
h(z) 0_ 



The entries of the polynomial matrix M(z) G C SX ( N ^ are of the form -kz r or constants; 
denoting A = z r , the matrices Mj(z), i = 1,2, can be expressed as 

Mi (A) = Mn - XM i2 

donde M u G C^-N+i)y-(N-i) y M<2i e c ( s -r+JV-i)x(JV-i)_ Ag a consequence) we have the 
following result: 
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Lemma 3 Assume that lankQ = r — N + 1. Then, rankG(z) = r for all z G C \ {0} i/ 
and on/y t/ rankM 2 (A) = N - 1 for all A G C \ {0}. 



Next step is to characterize when the rank of the matrix M 2 \ — XM 22 equals N — 1 for 
any A € C \ {0}. To this end, we use the Kronecker canonical form (KCF herafter) of the 
matrix pencil M 2 (A) (see the Appendix 6.2 for the details). By using the block structure 
notation A © B := diag( J 4, B), consider the KCF of the matrix pencil M2(A), i.e., 



K(A) := ^f*(A) © J M2 (A) © iV M2 (A) © S^f(A) 

where S^ ht (X) denotes the right singular part of M2(A), S^'(X) denotes the left singular 
part, Jm 2 (A) is the block associated with the finite eigenvalues of the pencil and, finally, 
Nm 2 (A) is the block associated with the infinite eigenvalue. Having in mind the structure 
of the different blocks appearing in the KCF of the matrix pencil M2(A), we can derive 
that the rank of K(A), and consequently of M2(A), is N — 1 for all A € C \ {0} if and only 
if K(A) has not right singular part and the only possibly finite eigenvalue is the zero one. 
In fact, we have the following result: 

Lemma 4 The rank of matrix M 2 (A) is N — 1 for each A G C \ {0} if and only if the 
following conditions hold: 

1. The KCF of the matrix pencil M 2 (A) has not right singular part and, 

2. If fj, is a finite eigenvalue of the matrix pencil M 2 (A), then fx = 0. 

Now, Lemma 4 allows us to decide when the rank of our initial polynomial matrix 
G(z) is r for all z G C \ {0}. Let us to remind all the given steps in reducing the initial 
polynomial matrix G(z): 



left. 



G(z) ~* G(z) ~» G(z) 



'Mi(z) G' 
M 2 (z) 



where 

1. G(z) = G(z)U(z), 

2. G(z)n r = G(z), 

3. G{z) = G(z)Q(z) 



l(z)\g], where Q G C 8 *<r- N+r > and Q(z) = diag[l, , 



l-r] 



4. If rank^ = r - N + 1, there exists R G C sxs invertible such that RG(z) 



where the matrix Q' G c( r - Ar+1 ) x ( r - iV + 1 ) is invertible, 

(A) = M a - \M i2 with A = z r 



h(z) 

5. The matrices Mj(z), i = 1,2, can be expressed as 
As a consequence, we have proved the following result: 



Theorem 2 Assume that supp&p C [0, N], where N G N with N > I, and take N < r < 
s. Let G(z) be the corresponding sxr polynomial matrix given in (6). Then, rankG(z) = r 
for any z G C \ {0} if and only if the following statements hold: 
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1. rankC? = r — N + 1 and, 



2. the KCF of the matrix pencil M2(A) has not right singular part, and the only possible 
finite eigenvalue is fj, = 0. 

For practical purposes it is not necessary to compute the KCF of the matrix pencil 
M2(A) (if possible). The needed information about M2(A) is obtained from its GUPTRI 
form (Generalized UPer TRIangular form). See the Appendix 6.3 for the details. As the 
matrix G(z) depends on z r , in what follows we identify the matrix G(z) with G(A) where 
A = z r . 



3.1 A toy model involving the quadratic B-spline 

The following example illustrates the result given in Theorem 2. Consider as generator (p 
the quadratic B-spline Ns(x), i.e., 

x 2 3 1 

N 3(x) = —X[0,1) + (-« + 3X - X 2 )X[1,2) + o( 3 - x fX[2,3) , 



where X[a,b) denotes the characteristic function of the interval [a, b). In this case, for the 
identity system Cf = f for all / G V v we have supp^C^ C [0,3], i.e., N = 3. Taking the 
sampling period T = 4/5, i.e., r = 4 and s = 5, the Laurent polynomials gi{z) given by 
(11) read: 



1 



1 



— z H — z 
2 2 

9 



37 



gs(*) 
»>(«) 

Following the above steps we obtain 



— z H 1 z 

50 50 25 

1 o 33 o 8 

— z ° H z 1 H z" 

50 50 25 



8 33 1 o 
g2(z) = 25 + 50 Z+ 50 Z 

2 -2 37 _! 9 
g4W = 25 Z + 50 Z +50 



G(z) 



r^ 4 

2^ 


±Z 5 

2^ 








33 4 
50^ 


1 y 5 
50^ 





8 r 3 
25^ 


25^ 





9 y 2 
50^ 


37 r 3 
50^ 





25^ 


37 r 2 
50^ 


9 .3 
50 ^ 


1 

- 50 


33 
50^ 


A z 2 

25^ 






Right multiplication by the matrix diagfl, z ,z 



1 .-2 z -3 } 



gives: 



G(A) = [M(A) | G] 



r ^a 


iA 


50 A 


i \ 

50 A 


2 X 
25 A 








2 

25 


1 


33 


- 50 


50 





± 

u 25 

JL 37 

50 50 

37 _9_ 

50 50 

^ 

25 u 
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where A 



The matrix Q £ C 5x2 has rank 2; performing some elementary operations 



on the rows of Q we obtain 





r 9 

50 


371 

50 




"0 





1 





0" 




"0 


0" 







8 
25 







1 
















8 

25 


G' = 










1 
















9 

50 


37 
50 












37 
9 


161 
18 





1 







37 
50 


9 

50 







0_ 




16 
L 9 


37 
9 








1 




8 

.25 


0_ 



= RQ 



Therefore, RG(X) = [RM(X) RQ] 



Ii(A) Q' 
I 2 (A) 



where 



\X 

5017 \ 
900 A 
_L _i_ H57 \ 
50 " r 450 A 



2 , 161 \ 
25 " r 900 A 



33 I 61 \ 
50 " r 450 . 



37 



In this case, a direct computation gives Km 2 (A) = Lj (A). As a consequence, Theorem 2 
ensures that the corresponding polynomial matrix G(z) possess a polynomial left inverse. 

Next we deal with the problem of computing a polynomial left inverse of G(z) whenever 
it does exist. 



4 Computing a polynomial left inverse of the matrix G(z) 

First note that if we compute a polynomial left inverse of the matrix G(A) then we obtain 
a polynomial left inverse of the matrix G(z). Indeed, remind that 

G(z) = G(z)U(«)fi- 1 Q(z) 

where U(z) = diagfz 7 *" 1 , (Wz) 7 "^ 1 , . . . , (W /r ~ 1 z) r ~ 1 ], Q r is the Fourier matrix of order r, 
and Q(z) = diagfl, z^ 1 , . . . , z 1_r ]. Thus, if L(z) is a polynomial left inverse of the matrix 
G(z), then the matrix 

L G (z) = diag[l, z,..., z r - l ]tt r diagfz 1 ^, (Wz) 1 ^, (W^ 1 zf-^z) 

will be a polynomial left inverse of the matrix G(z). Hence, we concentrate ourselves in 
computing a polynomial left inverse of the matrix G(z). To this end, consider G(A) = 
A T — XB T (A = z r ); being L(A) a polynomial left inverse of the matrix G(A), we have 
(.4 — Ai3)L T (A) = I r . Let us denote L(A) := L T (A). As we are searching for s x r matrices 
L(A), whose entries are polynomials, such that (A — XB)L(X) = I r we can use the following 
notation: 

L(A) = [Li(A) L2(A) . . . L r (A)] , i.e., Lj(A) denotes the ith column of L(A), 
Lj(A) = l° i +l}X + ---+t i X u , i = 1,2,..., r, where i\ € C s , k = 0,1,..., i/. 

As a consequence, (A — A,B)L(A) = I r is equivalent to 

A£°i + {At} - B^)X + ■■■ + (Al\ - Bt v ~ 1 )X v - Bl»X u+1 =P r , i = 1, 2, . . . , r, (13) 
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where P r denotes the ith column of the identity matrix I r . Equating coefficients, for each 
i = 1,2, . . . ,r, we obtain the set of linear equations 



Al\ = I\, Al\-Bl\ 



M\-Bt t 



0, -B£» = 0, 



or in matrix form 



—B 

A -B 

A -B 



A -B 
A 



U-l 



r 



, i = l,2,...,r, 



(14) 



where the resulting block matrix has order (v+2)r x (v+l)s. The goal is to find i/gN such 
that the above r linear systems become consistent. Next, we come back to the example in 
Section 3.1. 

The example revisited: Consider again the example involving the quadratic B-spline 
given in Section 3.1. In this case, G(z) = G(z)\J(z)£l^ 1 diag[l, z~ x , z 
A = z 4 we have 



2 , z 3 ] and, taking 



where 



A 



G(A) 



r^A 


\x 





33 \ 
50 A 


50 A 





25 A 





9 

50 





2 


37 


25 


50 


1 


33 


8 


- 50 


50 


25 





_8_ 

25 
37 
50 

JL 

50 



A T - XB T 







_8_ 

25 







_9_ 

50 
37 
50 



^ 

_2_ 

25 
37 
50 
_9_ 
50 



y B 



o 





_M 
50 
_ J_ 

50 







I 

"25 




































Here, the matrix S = 



-B 

A -B 
A 



of size 12 x 10 has rank 10. Choosing the columns of L(A) 



as Li (A) = £° + i]\ e C 5xl , the linear systems 



'-B 






\£ r 




"0" 


A 


-B 

A 




*■» 
£° 






. J 4. 



, z = l,2,3,4 



(15) 



have a unique solution. Observe that deleting the trivial equations 3 and 4, we have 
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consistent square systems. By using Matlab we obtain the left inverse 



4.4812 


-0.1438 


0.0166 


-0.0043 


-3.4840 


0.1118 


-0.0128 


0.0031 


1.6069 


-0.0514 


0.0056 


0.0000 


-0.4125 


0.0125 


0.0000 


-0.0000 


0.0500 


0.0000 


-0.0000 


0.0000 





"-0.0021 


0.0001 


-0.0000 


0.0000 




0.0517 


-0.0017 


0.0002 


-0.0000 


+ 10 3 


-0.4133 


0.0133 


-0.0015 


0.0004 




1.6071 


-0.0516 


0.0059 


-0.0015 




-3.4841 


0.1118 


-0.0129 


0.0033 



At this point, the challenge problem is to give conditions on the matrix pencil A — XB 
in order to obtain a left inverse with polynomial entries (having nonegative powers) by 
solving the corresponding linear systems (15). The answer to this question is based on 
the KCF of the matrix pencil A T — XB T . In our example the corresponding KCF is 
©iVi(A) (& L,J (A), i.e., the pencil has not finite eigenvalues, all the blocks associated 
with the infinite eigenvalue have order 1, and the left singular part has a unique block. In 
what follows, we prove that these conditions for the KCF of the matrix pencil G(A) are 
sufficient to give a positive answer to the raised problem in a very important particular 
case: 



4.1 The case where the oversampling rate is minimum for a fixed r > N 



It corresponds to the case where N < r and s - 
T = r/(r + 1). Here, the matrix pencil G(A) = A T 



r + 1, i.e., the sampling period is 
XB T has the form 



•• 


• 













0" 






* 


* 


• 


• 0" 


•• 


• 













* 






* 


* 


• 


• 


•• 


• 












* 




* 


* 





• 


• 


•• 


• 





* 




* 


* 


- A 


* 








• 


• 


•• 


• 


* 


* 




* 


* 




•• 


• 





• 


• 


•• 




* 


* 




* 







•• 


• 





• 


• 


* 




* 


* 







0. 




.0 •• 


• 





• 


• o. 



(16) 



i.e., denoting the entries of A T and B T by Ajj and Bjj respectively, we have Ajj = if 
2 + r<i + j<r + N + l, Bj N = and Bjj = if i + j > N + 1. Having in mind 
the structure of the matrices A T and B T we have rank(^4 T ) < r, rank(£> T ) < N — 1 and 
rank( [ _ s ] ) < r + N — 1. Whenever these matrices have maximum rank, the following 
result holds: 

Theorem 3 Assume that the singular matrix pencil A T — XB T of size (r + 1) x r satisfies 
the following conditions: 
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1. The pencil has not finite eigenvalues, 

2. rank(„4 T ) = r, 

3. rank(,B T ) =N -1, with N < r, and 
I rank([^ _ B ]) = r + N - 1 . 

Then, the Nr x (N - l)(r + 1) matrix 

-B 

A -B 
A 



-B 



A -B 
A 



has rank (N — l)(r + 1). 



First note that rank(>i T ) = r implies that the KCF of the matrix pencil A T — XB T 
has not right singular part (and also that is not an eigenvalue). Thus, by using Theorem 
2, the pencil A T — XB T has a polynomial left inverse. Before to prove Theorem 3, and 
in order to ease its proof, we first obtain, under the theorem hypotheses, the KCF of the 
matrix pencil A T — XB T : 



Lemma 5 The KCF of the matrix pencil A T — XB 1 



is 



© tfi(A) 

V 1=1 J 



J N- 



-i(A). 



Proof of Lemma 5: Since the matrix pencil has neither finite eigenvalues nor right 
singular part, we conclude that its KCF has the form iV(A) ©L Ze ^*(A), where N(X) denotes 
the blocks associated with the infinite eigenvalue and L Ze ^*(A) denotes the left singular part. 
Since r + 1 is the number of rows of the matrix pencil, r the number of columns, and the 
rank of B is N — 1 it cannot appear blocks of the form Lj (A) for i > N. Each left singular 
block increases in one the number of rows with respect to the number of columns; hence, 
as the size of A T — XB T is (r + 1) x r, it can appear only one left singular block in its 
KCF. Furthermore, we prove that this only left singular block corresponds to L^-_ 1 (A). 



Indeed, let /C^ 



A/Cg be the KCF of the matrix pencil A 



that rank(^l ) = rank(/Cj l ) = r, rank(,B ) = rank(/C^ 



TV 



rank 



The rank of the matrix 



-B 




A 


-B 


— Kb 






—Kb 



rank 



-Kb 



- XB . Obviously, we have 
1 and 

r + JV-1. 



coincides with its number of nonzero rows because the 



number of null rows of fCg is i — iV+ 1, i.e., the number of blocks in N(X); the matrix /C_4 

is 2r- (r-N + 1) = 



-Kb 



has not null rows so that, the number of nonzero rows of 
r + N-1. 

Assume that in the KCF of the matrix pencil A T — XB T appears a singular block 
Lj (X) with i < N — 1. Since the rank of B T is N — 1, the regular part in the KCF has a 
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block of the form iVj(A) with I > 2. By rearranging the blocks, we obtain that the KCF 
of A T - XB T is JVj(A) • • • © Lj(X); therefore 



-Kb 

K-A —K-B 




-1 

* * 



* * 
1 
1 

* * 








* * 



* * 











-1 

* * 















* * 



* * 



In this case, the rank of 



-Kb 

K-a —Kb 



is strictly smaller than r + N — 1 because the second 

row and the (r + l)th row are linearly dependent. This contradicts the hypotheses and, 
hence, the only left singular block is Lj f _ 1 (X). Having in mind that rank(£> T ) = N — 1, 

/r-N+l \ 

we conclude that the KCF of the matrix pencil A T - XB T is (J) N±(X) J © L^ r _ 1 (A). 

Proof of Theorem 3: Once we have determined the KCF of the matrix pencil 
A T — XB T we compute the rank of the matrix G r . If 1C_a — A/Cg is the KCF of the matrix 
pencil A — XB, it is obvious that 



rank(G r ) = rank 



-Kb 

K-A -K-B 

K-A -K-B 



K-A -K-B 
K-A. 



As /C^ — XKq is the KCF of the matrix pencil A T — XB T , Lemma 5 gives 



K 



L\ 



ICr. 




Ll 



where / = I( r _jv+i) denotes the identity matrix of order r — N + 1, and 






• 


•• 0" 




'1 


• 


•• 0" 


1 


• 


•• 


gC atx(JV-i ); l t = 





1 • 


•• 





• 


•• 







• 


• • 1 





• 


• • 1 







• 


•• 



£Nx(N-l) 
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As a consequence, 



rank(G r ) = rank 
















-L B 






I 














La 





-L B 






I 












La -L B 




L A 
I 






-L B 


L A 



A suitable interchange of rows and columns gives 



rank(G T 



rank 











-L B 
La 





-L B 



La 




-L B 
La. 



where the first r — N + 1 = r — rank(/Cg) are null rows; hence, the rank of G r equals 
(N — l)(r + 1) if and only if the remaining (N — l)(r + 1) rows are linearly independent. 
This is equivalent to the matrix 



£>A,B = 



-L B 
La 





-L B 



La 




-L B 
La 



^N(N-l)xN(N-l) 



has full rank. To prove it, we use the following result in [11, p. 32]: Let x(A) be a nonzero 
vector having the form 



x(A) = x + Axi + A 2 x 2 H h A £ x e , Xj e 



iNxl 



such that (La — ALg)x(A) = 0. Then, necessarily, e > N — 1. Now, let us continue by 
contradiction, and assume that the matrix Ca,b has not full rank. Then, there exists a 
nonzero vector z € C^^ -1 ^ 1 such that £a,b z = 0. Denoting z T = [zj [ _ 2 ■■■ zj Zq] 
where Zj G C Nxl , we obtain that 

(L A - XL B )(z + Azi + Az 2 + • • • + A^ 2 zjv_ 2 ) = 0, 
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which contradicts the minimal property for N — 1. Therefore, the matrix £-a,b has full 
rank and, finally, rankG r = (N — l)(r + 1). □ 

Remark: Note that Theorem 3 remains valid for any singular matrix pencil A T — XB T 
of size (r + 1) x r substituting N — 1 by p € N which satisfies < p < r. 

Consider the matrix pencil G(A) = A T — XB T of size (r + 1) x r with N < r. Assuming 
that the G(A) has polynomial left inverses, the following result gives sufficient conditions 
for computing one of such polynomial left inverses. Once we have got one solution, it is 
straightforward to derive the remaining solutions. 

Corollary 1 (Computing a polynomial left inverse of G(A)) LetG(X) = A T — XB T 
be a singular matrix pencil of size (r + 1) x r with N < r. Assume that G(A) admits 
polynomial left inverses, and that the following conditions hold: 

1. rank(„4 T ) = r, 

2. rank(fi T ) = N - 1, and 

3. rank([^ _ g ]) = r + N - 1 
Then, the linear systems 



~-B 






"0" 


A -B 




rp N ~ 2 i 






A -B 




fN-3 







A -B 




LS J 













p 



i= 1,2,... ,r, 



(17) 



where P r denotes the ith column of the identity matrix I r , admit a unique solution. More- 
over, let [i^ 2 if' 3 . . . £°] T € C^- 1 )^ 1 ) be this solution for i = 1, 2, . . . , r, and consider 

the polynomial vector Lj(A) = i°i+i}M Vtf~ 2 \ N ~ 2 , i = 1,2, . . . ,r. Then, the (r+l)xr 

polynomial matrix 

L(A):=[L 1 (A)L 2 (A)...L r (A)] 



satisfies 



L (A)G(A) 



Proof: Theorem 3 implies that the rank of the coefficient matrix G r £ (£ Nrx ( N - 1 )( r + 1 ) 
is (N — l)(r + 1) in (17). Having in mind (16), the last r — + 1 rows of B are null. 
Deleting these rows in the first row block (which become trivial equations in (17)), we 
obtain an square invertible matrix, and consequently (17) has a unique solution for each 
i = 1,2, ... ,r. Recalling (14), we finally obtain that L T (A) is a polynomial left inverse of 
G(A). □ 

Observe that any other polynomial left inverse A(X) of the matrix G(A) is given by 
A(X) = L T (A) + B(X) [I r+1 - G(A)L T (A)] , 
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where B(X) is an arbitrary r x (r + 1) polynomial matrix. 

For the matrix pencil G(A) = A T — XB T of size (r + 1) x r with N < r, it is easy to 
give sufficient conditions in order to satisfy the conditions 1-3 in Corollary 1. Namely: 

Corollary 2 Consider the singular matrix pencil G(A) = A T — XB T of size (r + 1) x r 
with N < r. Denoting A T = [Ajj] and B T = \Bj], assume that the following conditions 
hold: 

r + N + l 



A-j ^ if i + j = r + 2 or i + j 
Bjj^O ifi + j = N + landi>2, 

Then the conditions 1-3 in Corollary 1 are satisfied. 



(18) 
(19) 



Proof: Conditions (18) and (19) say that the entries marked as • in the matrices below 
are nonzero 



A = 









• 








• * 
• * * 



• 
• * 



* • 

• 







A u A 12 
A21 A22 



-<rx(r+l) 





* 


* 




• 


• 


• ' 






* 


• 





• 


• 


B = 


* 


• 


• 





• 


• 







• 


• 





• 


• 







• 


• 





• 


• 



B n 




^rx(r+l) 



where A 22 £ c( r - JV+1 ) x ( r - iV + 1 ) and B u G c^ N -^ xN . Trivially, rank(^4 T ) = r and 
rank(,B T ) = N — 1. Condition 3 comes by observing the form of the matrix [ -#]■ 
Interchanging rows and columns we obtain that the matrix [ _ B ] has the same rank 
than the matrix 

" 0" 
Bn 
A n B n A 12 
A 2 i A 22 
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Since the matrix A22 € C*- r N + 1 ) x ( r N + 1 ) [ s invertible, elementary row operations give 
the new matrix 















Bn 











An 


Bn 








A21 





-422 






Finally, the above matrix has rank 2(N — 1) + r — N+l=r+N— 1. 

□ 

Remark that condition 1 in Theorem 3 can be checked by using the algorithm guptri. In 
case that conditions 1-3 in Corollary 1 are satisfied, we could check directly the consistency 
of the linear systems (17); if they are not consistent, we derive that the pencil G(A) has 
not polynomial left inverses. 

5 Conclusion 

Consider the problem of the recovery of any function / in a shift-invariant space V v from 
the sequence of samples {(£/)(rn/ 's)} ne i of its filtered version Cf, where the positive 
integers r and s satisfy s > r, i.e., the oversampling setting. The existence of compactly 
supported reconstruction functions for this sampling problem is intimately related to the 
existence of a polynomial left inverse for a polynomial matrix G(z) associated with the 
sampling problem. This is equivalent to that the matrix G(z) has full rank for any z € 
C \ {0}. Other characterizations can be found in the mathematical literature involving 
the Smith canonical form of the polynomial matrix G(z) or the Euclides algorithm for 
the minors of order r in G(z). Unfortunately, whenever the parameter r is large, the 
aforesaid methods are useless from a practical point of view. In this work, by assuming 
that N < r < s, where dp C [0, N], we derive a new characterization for the existence of 
polynomial left inverses which involves the Kronecker canonical form of a singular matrix 
pencil. The advantage, from a practical point of view, of this new method is that we 
can retrieve the needed information from the KCF by using the so-called guptri algorithm. 
Furthermore, in the important case where s = r + 1, i.e., the oversampling rate is minimum 
(for a fixed r), we propose a method for the computation of a polynomial left inverse (and 
hence, the whole set of polynomial left inverses) of G(z). 

6 Appendix 

For the sake of completeness we include here a brief reminder on the canonical forms 
alluded throughout the paper. 

6.1 Smith canonical form of a polynomial matrix 

Recall that any mx n (m > n) polynomial matrix H(z) with rankH(z) = r (recall that 
the rank of a polynomial matrix is the order of its largest minor that is not equal to 
the zero polynomial) can be written as the product H(z) = V(z)S(z)W(z) where *V(z) 
an W(z) are unimodular matrices (i.e., the determinants of V(z) and W(z) are nonzero 
constants) of size m x m and n x n respectively, and S(z) is a diagonal mx n polynomial 
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matrix S(z) := diag[ii(z), . . . , i r (z), 0, . . . , 0]. Moreover, the diagonal entries (the so-called 
invariant polynomials of M(z)) are given by ij(z) = dj{z) / dj-\{z), j = 1,2, ... ,r, where 
dj(z) is the greatest common divisor of all minors of H(z) of order j, j = 1,2, ... ,r, and 
do(z) = 1. The matrix S(z) is called the Smith canonical form of the matrix H(z). See 
[15] for the details. 

6.2 Kronecker canonical form of a matrix pencil 

The Kronecker canonical form (KCF) for matrix pencils H(A) = A — \B, A, B G c mxn , is 
a generalization of the Jordan canonical form to matrix pencils (see [11]): There exist two 
nonsingular matrices U G £ mxm an d y g £nxn suc ] 1 ^ n ^Iq^ structure notation): 



— 1 qright 



(A)e j H (A)e7v H (A)e4 e/t (A), 



where Jh(A) © A r e(A) is the regular part of the matrix pencil, S^ ght {\) is its right singular 
part, and S l ^ t (X) its left singular part. The block Jh(A) is associated with the finite 
eigenvalues of the matrix pencil, and it reads: 



«/h(A) = Jhim) 



J l 9q (Mg 



where Jij(ni) is a Zj x Jordan block associate with the finite eigenvalue fii, i.e., 



1 



1 

Mi. 



1 



Recall that ^ is a finite eigenvalue of the matrix pencil H(A) if rankH(^) < rankH(A), 
being rankH(A) the order of the largest minor that is not equal to the zero polynomial. 
The block iVjgi(A) is associated with the infinite eigenvalue (if does exist) and it reads: 

A%(A) = iV Pl (A)©---©AT psoo (A) 

where N Pi G C PlX ' Pi is the matrix 



N Pi (A) 



and denotes the geometric multiplicity of the infinite eigenvalue which corresponds to 
the number of Jordan blocks for the infinite eigenvalue. Recall that the pencil H(A) has 
the infinite eigenvalue if its dual pencil H"(A) := AH(1/A) has the zero eigenvalue. 

If m ^ n or det(.A — XB) = for all A G C, then the matrix pencil also includes a 
singular part, S^ ght (X) and/or (A), and we say that the matrix pencil is singular. For 
the right singular part, we have 



aright 



(A) = L £1 (A) 
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where L £i (X) is a block of size x (gj + 1) defined by 



"0 1 



"1 



(20) 



1 



1 



Lq is a block of size 0x1 which contributes to a column of zeros. Analogously, the left 
singular part has the form 



where Li. (A) is a block of size + 1) x and Lq" is a block of size 1x0 which contributes 
to a row of zeros. 



6.3 GUPTRI form 



The GUPTRI form (Generalized UPer TRIangular form) for singular matrix pencils was 
done by Van Dooren [19, 20] by using unitary equivalence transformations. It is a gener- 
alization of the Schur-staircase form for matrices. 

Given a singular matrix pencil H(A) = A — AH with A, B G C mx ", there exist two 
unitary matrices U and V of size m x m and n x n respectively such that 



t/(4 - \B)V H 



A r — XB r 
A 




AH 



reg 







4 - AHj 



where the rectangular block upper triangular A r — A.B r and A\ — AH/ give the right and 
left singular structures of the matrix pencil, respectively. The remaining square upper 
triangular A reg — XB reg contains all the finite and infinite eigenvalues of H(A). Furthermore, 
the regular part A reg 



XB reg is in staircase form: 



A 



Aq * * 
Af -k 

An 



Breg — 



Ho * * 




where Aq — XBq and — AHoo reveal the Jordan structures of the zero and infinite 
eigenvalues, and Af — AHj, in generalized Schur form, includes the finite but nonzero 
eigenvalues. 

The eigenvalues fii are computed as pairs of values, denoted by (aj,/3j), ctj in the 
diagonal of A reg and 0i in the diagonal of B reg as follows: If a j 7^ and 7^ then 
/Xj is the finite nonzero eigenvalue jii = a.i/fii\ if on = and /3j 7^ 0, /Xj is the zero 
eigenvalue and; if «j 7^ and $ = then /Xj is the infinite eigenvalue. The case ai = 
and = does not correspond to an eigenvalue, instead it belongs to the singular 
part of the matrix pencil. In [8, 9] is described an efficient algorithm for computing the 
GUPTRI form of a matrix pencil. The implementation of this algorithm can be found in 
http : //www. cs .umu. se/research/nla/singular_pairs/guptri. 
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